Draft version October 14, 2008 

Preprint typeset using I^T^X style cmulateapj v. 10/09/06 



oo 
O 
O 
(N 
•*-> 
O 

O 



6 



> 

in 
oo 
(N 

o 

oo 
o 



X 



THE X-RAY HALO OF CEN X-3 
Thomas W. J. Thompson 1 & Richard E. Rothschild 1 

Draft version October 14, 2008 

ABSTRACT 

Using two Chandra observations we have derived estimates of the dust distribution and distance 
to the eclipsing high mass X-ray binary (HMXB) Cen X-3 using the energy-resolved dust-scattered 
X-ray halo. By comparing the observed X-ray halos in 200 eV bands from 2-5 keV to the halo profiles 
predicted by the Weingartner & Draine interstellar grain model, we find that the vast majority («70%) 
of the dust along the line of sight to the system is located within about 300 pc of the Sun, although the 
halo measurements are insensitive to dust very close to the source. One of the Chandra observations 
occurred during an egress from eclipse as the pulsar emerged from behind the mass-donating primary. 
By comparing model halo light curves during this transition to the halo measurements, a source 
distance of 5.7 ± 1.5 kpc (68% confidence level) is estimated, although we find this result depends 
on the distribution of dust on very small scales. Nevertheless, this value is marginally inconsistent 
with the commonly accepted distance to Cen X-3 of 8 kpc. We also find that the energy scaling of 
the scattering optical depth predicted by the Weingartner & Draine interstellar grain model does not 
accurately represent the results determined by X-ray halo studies of Cen X-3. Relative to the model, 
there appears to be less scattering at low energies or more scattering at high energies in Cen X-3. 

Subject headings: X-rays: ISM — pulsars: individual (Cen X-3) 



1. INTRODUCTION 

X-ray halos, which appear as diffuse emission sur- 
rounding X-ray sources, are created by small-angle scat- 
tering of soft X-ray photons from dust grains in the in- 
terstellar medium (ISM). Their study can provide infor- 
mation on interstellar grain properties (density, morphol- 
ogy, composition) and on the spatial distribution along 
the line of sight. Given a dust distribution, variability in 
the X-ray halo can be used to geometrically derive esti- 
mates of the distances to X-ray sources based on the time 
delays of the photons scattere d along the line of sight 
(jTriimper fc Sc honfclderl Il973h . Th is method has been 
applied in a number of cases fe.g.. iPredehl et all [2 000: 
Thompson et al.l 120061 : lAudlev et all 120061 : Kiang et all 
20071) . The goal of this work is to estimate the dust dis- 
tribution and distance to the eclipsing HMXB Cen X-3 
using measurements of its dust-scattered halo with the 
Chandra X-Ray Observatory. 

The theory of X-ray scattering and the production of 
X-ray halos has been described by a number of authors 
(e.g., Mathis & Lee 1991, Smith & Dwek 1998, Draine 
2003). Various interstellar grain models have been pro- 
posed. Mathis et al. (1977) developed a model com- 
posed of silicate and graphite grains with a size distri- 
bution of n(a) oc a~ 3 ' 5 , which reproduces the observed 
extinction of starlight. Weingartner & Drain (2001, here- 
after WD01) produced a grain model that additionally 
accounts for the diffuse infrared and microwave emission 
emission from the ISM by including sufficient small car- 
bonaceo us grains and polyc yclic aromatic hydrocarbons 
(PAHs). Zub ko et al.l (|2004T ) created a series of models by 
simultaneously fitting the extinction, infrared emission, 
and elemental abundance constraints, by including, for 
example, amorphous carbon particles, organic refractory 
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material, water ice, and voids. 

One of the many challenges in developing a viable 
model is that the characteristics of the dust may vary 
in different Galactic locations due to different local evo- 
lutionary histories, ISM phases, and abundances of ele- 
ments. Given such difficulties, a preeminent interstellar 
grain model has not been established. For our purposes, 
we use the WD01 grain model, although we allow the 
energy dependence of the scattering cross section to be 
a free parameter. 

The basic quantity that determines the characteristics 
of X-ray halos is the differential scattering cross section 
da/dn. This can be calculated using the exact Mie so- 
lution or the simpler Rayleigh-Gans approximation. For 
the WD01 grain model, the differential scattering cross 
section as a function of scattering angle 4> can be approx- 
imated by the simple analytic form 
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where </> s ,50 ~ 360" (keV/E) is the median scattering an- 
gle as a function energy (Draine 2003). It is important 
to note that the following results are probably somewhat 
dependent on this chosen form of the differential scatter- 
ing cross section. 

Cen X-3 is one of the brightest accreting X-ray pul- 
sars, and it is one of the six pulsars in which the ob- 
servation of eclipses has permitted the determination of 
all orbital and stellar parameters. The rotation period 
of the pulsar is 4.8 s, and the binary orbital period is 
2.1 days. The inclination of the orbital plane and the 
mass of the neutron star have been estimated to be 
i ~ 70° and M = 1.21 ± 0.21 M (Ash et al. 1999). 
The optical counterpart of Cen X-3 has been identified 
as an early t ype star of radius -R ~ 12 Rq and mass 
M ~ 20 M Q (Rappaport k. Josslfl983] ). Accretion in the 
system probably occurs via an accretion disk (Bildsten 
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1997; Takeshima et al. 1992), although the mass trans- 
fer mechanism probab ly includes an X-ray excited wind 
(pav fc St even g fl993h . 

Although well-studied, the distance to Cen X-3 re- 
mains uncertain. Using a BO II stellar classification 
model, an 8 kpc distance to Cen X-3 was obtained 
(|Krzeminskil 1 1 974} ) . however, subsequent work indicated 
that the supergian t primary was of type 6-8 III (Hutch- 
ings et al. 1979). IDav fc Tennanti (|1991[) observed soft 
emission during two eclipses using EXOSAT and at- 
tributed it to the dust-scattered halo, from which they 
derived a distance of 5.4 ±0.3 kpc. Discrepancies of ^2.5 
kpc in distance estimates to Cen X-3 lead to nearly a 
factor of two difference in the inferred optical and X-ray 
luminosities. 

Cen X-3 (I = 292? 09, b = 0?34) is particularly suited 
for X-ray halo studies. The source is usually very bright, 
and the interstellar hydrogen colum n density along the 
line of sight is 1.20 x 10 22 cm~ 2 (Dic key fc Lockmanl 
1990), implying a la rge enough optical depth to result in 
an appreciable halo (|Predehl fc Schmittlll995l ) , yet not so 
large as to include a substantial contribution due to mul- 
tiply scattered photons. Moreover, the eclipsing nature 
of the binary system is advantageous because it should 
produce the largest fraction of variability in the X-ray 
halo, which makes a geometric distance measurement 
more direct. 

Section 2 presents the general equations that char- 
acterize the shape and variability of X-ray halos. For 
our study of the X-ray halo of Cen X-3, two Chan- 
dra High Energy Transmission Grating (HETG) obser- 
vations are used; one observation took place outside of 
eclipse (<i> or b ~ 0.4 — 0.6) when there was little variabil- 
ity in the flux (hereafter, called the "plateau" phase of 
the orbit based on the shape of the light curve), and 
the other occurred during an egress from eclipse (§ [3]). 
Because accurate halo measurements critically depend 
on careful subtraction of the instrumental point spread 
function (PSF), the shape and normalization of the PSF 
is parametrized using all Chandra /HETG observations 
of the unabsorbed sources Her X-l and PKS 2155-304. 
These results are presented in the appendix to the pa- 
per. The method for deriving the dust distribution to 
Cen X-3 is discussed in § 14.11 With a dust distribution 
established, the large amplitude change in flux during 
eclipse egress is used to constrain the distance to the 
system (§ 14. 2p . A comparison of the inferred dust distri- 
bution to observations of CO emission, star counts, and 
interstellar extinction in the direction of Cen X-3 is pre- 
sented in § 14.31 We discuss the implications of our results 
in § 5, and we provide a brief summary in § [5] 

2. X-RAY HALOS 

Following Draine & Tan (2003; hereafter DT03), for a 
steady source, the intensity of single-scattered photons 
arriving at halo angle 9 is given by 

h{e)=F^T sc J 1 {9), (2) 

where Fx and r sca are the flux and scattering optical 
depth, respectively, where 

" p(x) da 



h(0) 



-dx. 



(3) 



angle are related through cf) w 0/(1 — x), and the nor- 
malization is chosen such that J I 1 {9)2Ti6d6 — 1. For a 
variable source, a photon appearing at halo angle 9 after 
scattering from a dust grain at a fractional distance x to 
a source D kpc away will be delayed with respect to the 
central source by 



5t = 1.21- 



(— y 

\ arcsec / 
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kpc V arcsec J 1 



(4) 



l (l-x) 2 dQ 

p(x) is the dimensionless dust density at fractional dis- 
tance x to the source, the scattering angle and the halo 



In this case, Fx must be moved inside the integral (eq. 

3) because the halo flux at time t is proportional the 
flux of the central point source at the "retarded time" 
£' = t- St, i.e., h{t,6) oc F x (t - 8t(x,6)). 

It is clear that to accurately determine the proportion- 
ality between the source flux and the halo, one must know 
the history of the source flux for a sufficiently long time 
period preceding the observation. In most cases this is 
not practical, so typically one assumes that any change 
in the source flux took place sufficiently long ago for the 
halo to have completely responded to it (as we do when 
estimating the dust distribution, § 14. ip . or one models 
the history of the source flux using reasonable estimates 
(as we do when constraining the source distance, S 14. 2|) . 

The intensity of the X-ray halo due to photons scat- 
tering two or more times (I2 , I3, ■ ■ ■ ) is straightforward 
for a steady source given an assumed dust distribution, 
and we refer the reader to DT03 for the appropriate re- 
cursion formula. However, accounting for the time delays 
for multiply-scattered photons is extremely cumbersome. 
Fortunately, the intensity of the halo due to photons 
scattered n-times, relative to the observed source flux, 
is T™ a /n! (jMathis fc Led 1 199 if) , so the halo intensity for 
multiply-scattered photons decreases rapidly for reason- 
ably small optical depths (r sca < 2). Usually, the first- 
and second-order scatterings are sufficient to model the 
total halo. 

3. OBSERVATIONS & DATA EXTRACTION 

In this paper, we use two Chandra/HETG observations 
of Cen X-3 (ObsIDs 1943 & 7511), and several Chan- 
dra/HETG observations of PKS 2155-304 and Her X-l 
to model the PSF (see the appendix). Although the use 
of the grating reduces the zeroth-order image count rate 
by about a factor of 3, the dispersed spectrum provides 
accurate measurements of the source flux and spectrum, 
which is necessary to accurately model the character- 
istics of the X-ray halo and the PSF. One of the ob- 
servations (ObsID 1943) began at MJD 51908.01 and 
lasted for 45 ks, corresponding to $ or b ~ 0.38 — 0.63, 
where determination of the orbital phases was accom- 
plished using the mid-eclipse ephemeris from Burderi 
et al. (2000), and the orbital period and evolution 
[Potb/Porb = -(1-738 ± 0.004) x 10~ 6 yr- 1 ] from Na- 
gase et al. (1992). A dip in the flux was observed at 
about <& or b ~ 0.57, so we ignored data after this time to 
simplify analysis by removing any time dependence (eq. 

4) . X-ray halo measurements from this observation were 
used to estimate the dust distribution along the line of 
sight. The other observation of Cen X-3 (ObsID 7511), 
which resulted from a successful observing proposal with 
the primary goal of geometrically estimating the dis- 
tance to the system, took place during an egress from 
eclipse. The 40 ks observation began at MJD 54355.46 
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(^orb ~ 0.08 — 0.30). The preceding mid-eclipse epoch 
was calculated to be MJD 54355.313 ± 0.007. Figure Q] 
(top-left panel, solid curve) shows the Chandra/HETG 
light curve during the eclipse egress observation. We 
also present an RXTE/PCA light curve for comparison. 
Clearly, the Chandra observation occurred during a rela- 
tively lower flux state, leading to a more gradual increase 
in flux coming out of eclipse as compa red the PCA obser- 
vation. iPriedhorskv fc Terrelll (|1983[ ) discussed the ape- 
riodic 120- 165 day time s cale b etween the different flux 
states, and IClark et alj (|1988l ) studied how the stellar 



wind and the shape of the egress light curves are affected 
by the different states. 

Data analysis was performed using the standard tools 
of the Chandra Interactive Analysis of Observations 
(CIAO) software version 4.0 and Calibration Database 
(CALDB) version 3.4.3. A single first-order Medium En- 
ergy Grating (MEG) and High Energy Grating (HEG) 
spectrum was extracted for the plateau phase observa- 
tion (ObsID 1943). For the eclipse egress observation 
(ObsID 7511), ten 3968 s first-order MEG/HEG spec- 
tra were extracted to model the evolution of the point 
source spectrum. Spectral models were fit to each spec- 
trum using XSPEC version 11.3, and fluxes (in units of 
[photons/cm 2 /s]) were measured in 15 energy bands of 
200 eV width from 2-5 keV. Finer time resolution on the 
flux evolution was obtained by comparing the average 
HETG count rates in each time range and energy band 
to the relative count rates in smaller 100 s time bins. 

Separate X-ray halo images and exposure maps were 
extracted in 200 eV bands from 2-5 keV, from which 
we obtained exposure-corrected surface brightness distri- 
butions using 40 logarithmically-spaced annular regions 
surrounding the source. For the plateau phase observa- 
tion, a single exposure map was created for each energy 
band. On the other hand, accurate measurements of the 
evolution of the halo flux in each energy band during 
the eclipse egress observation required the use of sepa- 
rate exposure maps (due to the dither of the telescope) 
using 500 s bins (1200 exposure maps in total). Due 
to issues with the diffracted HETG photons confusing 
the background analysis (HETG photons create a dif- 
fuse "transfer swath" analogous to the "transfer streak" 
in the zeroth-order image), we did not use halo angles 
greater than 110". The surface brightness distributions 
were divided by the flux measurements in each band to 
produce images in units of flux fraction per square arc- 
second. Halo angles smaller than about 3" are affected 
by pile-up and ignored. To minimize concerns due to 
multiply-scattered photons, we excluded energies below 
2 keV in each observation. 

In the appendix, we parametrize the PSF as the frac- 
tion of the point source flux comprising the PSF as a 
function of halo angle. Although the observed surface 
brightness distributions represent the convolution of the 
PSF and the X-ray halo, in our analysis we treat it as the 
sum. This is acceptable because we restrict our investi- 
gation to angular scales that are much larger than the 
~0'/5 resolution of the telescope mirrors. The PSF con- 
tribution in each energy band is simply subtracted from 
the observed surface brightness distribution. The result- 
ing X-ray halos, rebinned in 1 keV bands, are shown in 
Figure [2] The fitted halo components are discussed be- 
low. 



4. ANALYSIS & RESULTS 
4.1. Estimating the Dust Distribution 

In order to estimate the dust distribution to Cen X- 
3, we used the plateau phase Chandra observation. We 
assumed the halo had sufficient time to respond to the 
eclipse egress, which based on the orbital ephemeris oc- 
curred ~30 ks earlier. Once we established a dust distri- 
bution, we checked whether or not the assumption was 
acceptable. First, we created separate single-scattered 
X-ray halo profiles for twenty step-function dust distri- 
butions of width Ax = 0.05 (hereafter referred to as 
"clouds"), spanning the entire line of sight, i.e., 0.00 to 
0.05, 0.05 to 0.10, and so on, for each energy band. By 
summing the twenty profiles with equal weighting, one 
could obtain the single-scattered X-ray halo for dust dis- 
tributed uniformly along the line of sight. The normal- 
ization of each cloud (a x ), which is proportional to the 
contribution to the total scattering optical depth (due 
to the normalization J Ii(6)2ivQdO = 1), was allowed to 
vary freely. The initial fit function at each energy is 



%^ = f>lf(0), 



x=0 



where 



,0.06(*+l) 1 d(J 



(5) 
(6) 



Note that the superscript in If is a label that indicates 
a position along the line of sight and is not an exponent. 
Using this method, the dimensionless dust density and 
scattering optical depth are simply 



p(x) 



(7) 



respectively. The halos for each energy band were fit 
simultaneously. For the initial fit, we used the energy 
scaling of the scattering optical depth from the WD01 
interstellar grain model (i?y = 3.1); however, we found 
it to underpredict the halo flux at higher energies (or 
overpredict the flux at lower energies). We therefore al- 
lowed the energy scaling to also be a free parameter, 
although wc required that it be a smooth function, i.e., 
Tsca(E) oc [E — E )~P, where E and (3 are fit parame- 
ters. Nevertheless, the angular dependence of the scat- 
tering at a particular energy still follows the functional 
form predicted by the WD01 model (see eq. [1]). The 
best-fit normalization of each cloud provided a prelimi- 
nary dust distribution. Given the preliminary model, the 
X-ray halo for double-scattered photons was calculated 
(12(6)), included in the model as an additional term in 
eq. (5), and the normalization for each cloud was refit. 
This process was repeated in an iterative manner to cov- 
erage to the dust distribution whose single- and double- 
scattered X-ray halos accurately described the Cen X-3 
halo. Third- and higher-order scatterings were ignored, 
but their contribution to the halo is negligible. The Cen 
X-3 X-ray halos in units of source flux and the best-fit 
halo models are presented in Figure[2j and the integrated 
(2-5 keV) surface brightness distribution is shown in Fig- 
ure [31 Evidently, the 3"-100" halo can be modeled using 
only four separate dust clouds (Fig. [3J dotted curves). 

The energy dependence of the scattering optical depth 
is presented in Figure [4l showing the relatively brighter 
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Fig. 1. — Comparison of the eclipses and flux states in count rates (left panels) and fluxes (right panel) from 3-5 keV between the 
Chandra eclipse egress observation (ObsID 7511) and an RXTE/PCA observation. No data were o btai ned for the Chandra eclipse ingress 
(dotted curve). The ingress is modelled as the mirror image of the egress (multiplied by 0.87, see £j 14. 21 . 
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Fig. 2. — Cen X-3 halo measurements and best-fit model halos (§ 14.11 shown in three bands from 2-5 keV in units of source flux per square 
arcsecond. Note that the actual fits to the halos used 200 eV bands. Larger 1 keV bins are displayed to facilitate the visual comparison 
between the models and the data. The dotted curves show the contribution to the halo from the individual single-scattered (Ax = 0.05) 
clouds (see Fig. |3]for the cloud locations), the dashed curves show the halo contribution from photons scattered twice, and the solid curves 
show the sum of the single- and double-scattered halos (xi = 1.05, 502 dof). Third- and higher-order scatterings are nearly negligible and 
are not included in the model. 



halo at higher energies as compared to the WD01 inter- 
stellar grain model predictions (see discussion). The Cyg 
X-l and GX 13+1 points are the same data as presented 
in Fig. 11 of Draine (2003), who compared the opti- 
cal depth calculated from the WD01 model to observa- 
tions. Although Draine (2003) extended his comparison 
to other sources, only Cyg X-l and GX 13+1 span the 
2-5 keV energy range that we investigate for Cen X-3. 

One of the primary goals of this work is to geometri- 
cally constrain the distance to Cen X-3, therefore it is 
important to quantify the uncertainty in the dust distri- 
bution because this directly affects the size of the time 
delays of photons scattered along the line of sight. By 
fixing the double-scattered halo model to the one cor- 
responding to the best-fit dust distribution, the single- 
parameter 90% confidence uncertainties in the contribu- 
tion from each of the twenty clouds was calculated. In 
each case the lower limit on the cloud size was zero, i.e. 
no dust, simply because the clouds to either side of the 
cloud in question could produce a surface brightness dis- 
tribution with a similar shape. On the other hand, the 
upper limits to the size of each cloud are well-determined. 



Figure [5] shows the best-fit dust distribution towards Cen 
X-3 and the single-parameter upper limits for the size of 
each cloud. Three examples of different dust distribu- 
tions that are also consistent with X-ray halo data are 
presented. The majority of the dust («70%) in each case 
is located at small fractional distances to Cen X-3. Note, 
however, that we do not have reliable data for halo angles 
less than about 3", so we are mostly insensitive to dust 
that produces centrally-peaked halos, i.e., dust very near 
the source (x > 0.99). Therefore, it is more accurate to 
consider the scattering optical depth in Fig. [4] as a lower 
limit. 

It could be argued that lack of knowledge of the dust 
very near the source poses a problem for our distance es- 
timate. This is not the case: For an assumed distance of 8 
kpc, the ~0'/5 angular resolution of the Chandra mirrors 
corresponds to approximately 4000 AU, which, based on 
the orbital parameters of the binary system, is more than 
4 orders of magnitude larger than the semi-major axis of 
the orbit. Therefore, the point source flux already in- 
cludes all of the effects due to scattering in the vicinity 
of the binary system. 
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Fig. 3. — Surface brightness distribution (with PSF subtracted) 
from 2—5 keV. As in Fig. [2] the dotted curves show the contribution 
to the halo from the individual (Ai = 0.05) clouds and the dashed 
curves show the halo contribution from photons scattered twice. 
The location of the dust for each single-scattered halo component 
is labelled at left (also see Fig. QjJ, 
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Fig. 5. — Best-fit model dust distribution to Cen X-3 (solid 
histogram). The majority of the dust is located near the Sun 
(0.00 < x < 0.05), and the other primary dust locations are at 
x > 0.75. Although the best-fit distribution accurately models the 
data (xl = 1.05, 502 dof), the size of any individual dust cloud 
is not well-constrained (dotted histogram). The three small panels 
labeled a, b, and c (corresponding to the upper limits shown in 
the main panel) show examples of dust distributions that are also 
consistent with the data at the 90% confidence level. 
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Fig. 4. — Best-fit scattering optical depth to Cen X-3 as a func- 
tion of energy (solid curve). The empirical best-fit curve is de- 
scribed by T sca (£) oc (E - Eq)-! 3 , where E ~ 1.96 and /3 ~ 0.32. 
The optical depth predicted by the WD01 interstellar grain model 
(Rv = 3.1) is shown as a dash-dotted curve. The normalization for 
the WD01 model was chosen so that the areas under each curve 
are th e same. Also included are Cyg X-l d ata from lYao et al.l 
J2003T> (filled circles) and GX 13+1 data from lSmith et al.l t206W 
(triangles) (see discussion). The normalizations for the Cyg X-l 
and GX 13+1 data have been modified to more closely align the 
points with the e mpirical curve for Cen X-3. The data points from 
lYao et al. (2003) have been correc t ed a ssuming uniform dust, and 
the data points from Smith et al. (2002) have been corrected as- 
suming the dust density is proportional to distance. See Draine 
(2003) for further details on the correction method. 

Considering the best-fit dust distribution, it is clear 
that the assumption that the halo during the plateau 
phase observation had responded to the preceding eclipse 
egress is not valid for all halo angles. As pointed out 
above, an estimate of the time Cen X-3 has been out 
of eclipse at the beginning of the observation is about 
30 ks. The causally-connected regions of x vs. 9 phase 
space for delay times less than 30 ks are shown in Figure 
[6] for three source distances, indicating that the assump- 



10 100 
Halo Angle (arcsec) 

Fig. 6. — Phase space boundaries (x vs. 8) showing causally- 
connected regions for delay times less than 30 ks. The hatched 
regions of phase space have yet to respond to changes in the source 
flux 30 ks earlier (see eq. 4). For example, assuming a source 
distance of 4 kpc and dust at x = 0.75, for 6 > 46" the halo has 
not yet responded to changes in the source flux 30 ks earlier (see 
point and arrow). 

tion is invalid for large x and 9. Fortunately, the vast 
majority of the dust along the line of sight is located at 
x < 0.05 which presents no inconsistency, and the two 
clouds closest to the source (0.90 < x) dominate the halo 
at 9 < 10" (see Fig. [3]), where the halo angles are small 
enough to have responded to the previous eclipse egress. 
The only potential problem concerns the x — [0.75,0.80] 
dust cloud, which for 9 > 30" have not had sufficient time 
to respond to source flux variability. However, the dust 
cloud at x < 0.05 also begins to dominate for sufficiently 
large halo angles, so we conclude that the assumption is 
reasonably valid. 

4.2. Constraining the Source Distance 
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Fig. 7.— RXTE/ ASM light curve of Cen X-3 for the 60 days 
(corresponding to x = 0.98 for D = 8 kpc and 8 = 100") preceding 
the Chandra observation during eclipse egress. The long-term Cen 
X-3 light curve was modeled using the Chandra fluxes scaled by the 
ASM count rate. Orbital phase 0.0 corresponds to MJD 54355.31. 

Because our goal is to use the Cen X-3 halo during 
the eclipse egress to constrain the distance to the sys- 
tem, it is necessary to model the source flux history for 
a sufficiently long time period preceding the observation 
(§ E|). We chose to model the flux history by first as- 
suming the shape of the light curve during the eclipse 
ingress in each energy band is same as the shape dur- 
ing egress, but scaled to 87% of flux after egress (Fig. [TJ 
top- left panel, dotted curve) . It is widely known that the 
Cen X-3 flux decreases after 3> or b ~ 0.5, and the scaling 
factor represents the average 2-10 keV flux change that 
was observed duri ng two consecutive binary orbits by 
ISuchv et all ((20081. So-called "pre-eclipse dips" are also 
commonly seen in Cen X-3 light curves (Fig. [TJ lower- 
left panel), however, we did not include these features 
in our model. Second, for the 60 days prior to the ob- 
servation we modeled the flux history by multiplying the 
Chandra light curve (and modeled ingress) by the relative 
RXTE/ ASM count rate at the time (linearly interpolat- 
ing the ASM rates at intervening times). The ASM light 
curve for Cen X-3 is shown in Figure [7] 

Admittedly, the modeled flux history is rather approx- 
imate. Not only have we ignored the pre-eclipse dips, 
but the assumption that the shape of the Cen X-3 light 
curve is independent of the absolute flux (the relative 
ASM count rates span ~0.05x to 5x of the flux dur- 
ing Chandra observation; Fig. [7|) is incorrect (see Fig. [TJ 
right panel). Yet most of the flux at sufficiently large halo 
angles is due to dust located within 300 pc of the Sun. 
The corresponding delay times are small, and the precise 
history of the source light curve can be measured directly 
from the Chandra observation. Therefore, any inaccura- 
cies in the modeled flux history will have a very limited 
effect, so a crude model may be assumed sufficient. On 
the other hand, at smaller angles the halo can be at- 
tributed to dust at larger x, making any inaccuracies in 
the flux history model more problematic. For this reason, 
we ultimately chose not to investigate halo angles smaller 
than ~20". Separate 200 eV PSF-subtracted surface 
brightness distributions were extracted, using time bins 
of 1.5 ks. The energy bands were summed together, and 
halo light curves were created for four annular regions: 
[22'.'8,33'.'5], [33"5,49"2], [49'.'2, 72'.'2], & [72'.'2, 106'.'0]. 

Using the model flux history, we created model halo 
light curves for various energies, halo angles, and source 
distances, using the best-fit dust distribution as well as 
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Fig. 8. — Fit statistic (x 2 ) versus distance to Cen X-3. At the 
68% (90%) confidence level, the distance is 5.7±°\l (1°'?) k P c - The 
quality of the fit at the best-fit distance is x 2 /dof = 87.0/96 = 0.91. 
A more appropriate method for error determination is discussed 
below. 

the three example dust distributions that are also consis- 
tent with the X-ray halo data. The final 2-5 keV single- 
scattered model light curves for the four annular regions 
are then derived by summing the energy bands and inte- 
grating over the annular regions, viz. 

14 Q 

h (t, A9) = f; 1 [ 2 27T0/1 (t, 6)dQ, (8) 
i=o ^ a ~ l > Je i 

where the superscript i is a label that indicates the center 
of the 200 eV energy band at (2.1 + 0.2?) keV, and 6>i and 
62 are the boundaries of each annular region. 

Although the single-scattered halo light curves can 
be easily modeled, we must also account for multiply- 
scattered photons, as well as our ignorance of the specific 
form of the dust distribution (Fig. [5]). Calculating the 
time dependence of multiply-scattered photons is compu- 
tationally cumbersome (DT03). Fortunately, these pho- 
tons comprise ~10% of the total halo flux over the angles 
of interest (Fig. [3]), so a rather simple treatment should 
be sufficient. Moreover, because the dust is concentrated 
at small x, the response of the halo to multiply-scattered 
photons should occur almost as rapidly as the response 
to single-scattered photons. Thus, to roughly account for 
multiply-scattered photons, we added 10% (the approx- 
imate value of I2/I1 over the angles of interest in the 
plateau phase observation) to the normalization of the 
single-scattered model light curves. We also allowed the 
normalization to vary by 10% when obtaining the fits. 

To address the uncertainties in the dust distribution, 
we calculated model light curves for the best-fit distribu- 
tion and for the three example distributions shown in Fig. 
[5l Although the general characteristics of the light curves 
(shape and normalization) for the different dust distribu- 
tions were similar, subtle differences were present. The 
size of the differences in the light curves were about 7% 
± 3% (relative to the model light curve for the best-fit 
distribution) at the beginning of the observation (and 
for the most part, independent of distance and angle), 
and 4% ± 1% at the end. Clearly, these measured dif- 
ferences only apply to the selected set of four potential 
dust distributions. In each case, however, the majority 
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Time (ks) 

Fig. 9. — Eclipse egress light curves from 2-5 keV for four annular regions (labeled in each panel). The model light curves for D = 5.7 
kpc are shown as solid curves. The dotted curves show the HETG light curve, and the dashed curve shows what the predicted light curves 
would be if the prominent dust cloud was shifted from x = [0.00,0.05] to x = [0.20,0.25]. 
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Fig. 10. — Histogram of the best-fit distances to Cen X-3 from 
2500 synthetic data sets. The distribution is fitted with the sum of 
two Gaussians. The standard deviation of the distribution is 1.45 
kpc, and the 90% confidence level lower and upper limits, which 
are determined from the model of the distribution, are 2.1 kpc and 
8.3 kpc, respectively. 

of the dust is located relatively close to the Sun, so the 
differences in the light curves from this set should pro- 
vide a decent measure of the uncertainty in l\ (t, AO) due 
to the uncertainty in p(x). We chose to account for the 
uncertainty in the dust distribution by adding ~4%-7% 
(depending on differences in the model light curves as a 
function of time) to the intrinsic uncertainty of the halo 
light curves. 

Finally, the distance to Cen X-3 was obtained by fitting 
the model halo light curves for the four annular regions 
and for source distances between 3 kpc and 12 kpc in 
steps of 0.1 kpc to the observed light curves. For each 
distance, the x 2 statistic was used to determine the qual- 
ity of the fit. Our results indicate that the distance to 
Cen X-3 is 5.7±g ;| (±° 9 7 ) kpc with 68% (90%) confidence 
(Figure [5]). Note that we derive a more reasonable esti- 
mate of the distance uncertainty below. 

The comparison of the measured halo light curves to 
the model fits for a source distance of 5.7 kpc (solid 
curves) is shown in Figure [5J Also shown is the model 
light curves if the x = [0.00, 0.05] dust cloud were moved 
to x = [0.20,0.25]. The halo light curves are therefore 
clearly consistent with most of the dust having a local 
origin. 
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Fig. 11. — Similar to Fig. \E\ but decreasing the width of the 
x = [0.00, 0.05] dust cloud to Ax = 0.01 (while keeping the total 
amount of dust in the cloud constant by multiplying the normal- 
ization by 5) and trying each position within x = [0.00, 0.05]. 

Because the model light curves are not strictly linear 
in distance, A^ 2 may not be a legitimate estimator of 
the distance uncertainty (Fig. [5]). A more reasonable 
estimate of the uncertainty can be obtained using Monte 
Carlo simulations. We accomplished this by simulating 
2500 synthetic data sets, assuming the light curve mod- 
els for D = 5.7 kpc accurately describe physical reality, 
and assuming the measurement uncertainties are nor- 
mally distributed. For each realization the best-fit dis- 
tance is determined from the x 2 minimum. The resulting 
probability distribution of best-fit distances is presented 
in Figure [TOl indicating that the uncertainty in distance 
based on A\ 2 is indeed too small, and that a more rea- 
sonable distance estimate to Cen X-3 is 5.7 ± 1.5 (l 2 'g) 
kpc with 68% (90%) confidence. The 68% error limits 
were determined simply from the standard deviation of 
the distribution, and the 90% limits were determined by 
integrating over the model of the distribution, which is 
composed of the sum of two Gaussians. 

Another issue to address is the use of Ax — 0.05 bins 
for the dust clouds. The uncertainty in the position of 
the dust within the x = [0.00, 0.05] cloud can have a sig- 
nificant effect on the model halo light curves. To test the 
effect of various dust positions for x < 0.05 we calculated 
model light curves for the five positions with Ax = 0.01 
in x < 0.05, maintaining a constant amount of dust in the 
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Fig. 13. — Histogram of the number of stars within 1?25 of the 
line of sight to Cen X-3. Data are based on the Hipparcos Catalog. 
The fractional distance assumes a source distance of 5.7 kpc. 



Fig. 12. — CO emission observed in the direction of Cen X-3 
(Dame et al. 2001), with the two most prominent emission features 
labeled. 

cloud by increasing its normalization by a corresponding 
amount. For each of the five trial dust positions, the fit 
statistic was calculated at each distance (Fig. [TT]) . The 
results are somewhat troubling, because although dust 
positions at x < 0.04 do not greatly affect the implied 
distance, dust oriented at x = [0.04, 0.05] produce two 
local minima, with the best-fit distance becoming very 
small. 

4.3. Comparison of Dust Distribution to CO Emission, 
Star Counts, and Interstellar Reddening 

The dust distribution to Cen X-3 can be checked for 
consistency in a number of ways. Here we apply three 
methods. First, using the results of Dame et al. (2001), 
we can map the CO emission in the direction of Cen X-3, 
which presumably also traces the dust. The most promi- 
nent peak in the emission has a velocity of —9.8 km/s 
(Fig. [12]). The kinematic distance to the CO responsi- 
ble for this feature could be estimated from the rotation 
curve of the Galaxy, but because molecular clouds have 
a cloud-cloud dispersion of 4—5 km s , kinematic dis- 
tances for clouds with \v\ < 10 km s -1 are not reliable; 
all that can be said is that CO emission has a relatively 
local origin (T. M. Dame, priv. comm.). Second, par- 
allax measurements provide the spatial density of stars 
along the line of sight to Cen X-3. Using data from the 
Hipparcos Catalog, we created a histogram of the num- 
ber of stars within 1?25 of the line of sight (Figure [T3|) . 
Of course, the parallax data are subject to certain selec- 
tion effects. For one, stars at larger distances are less 
likely to be above the flux limit of the survey. On the 
other hand, a given solid angle corresponds to a larger 
volume of space at larger distances. Nevertheless, as- 
suming star-forming regions are cospatial with regions of 
higher dust density, the histogram of star counts supports 
that conclusion that most of the dust is located at rela- 
tively small distances from the Sun. Finally, interstellar 
reddeni ng data may be a ble to constrain the dust distri- 
bution. iMarshall et alJ ()2006l ) modeled the interstellar 
extinction distribution in the Galaxy by comparing the 
infrared colors measured by the Two Micron All Sky Sur- 
vey (2MASS; Cutri et al. 2003 ) to the colors fro m the 
Besanoon model of the Galaxy (|Robin et al.| [2003). The 
resulting data for the Cen X-3 sightline, which shows the 
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Fig. 14. — Interstellar extinction (Afc s and &Ak s kpc -1 ) versus 
distance in the direction of Cen X-3 (cen t ered ? 13 away from the 
sightline). Data are from Marshall ct al. (2006). 

interstellar extinction (Ak b ), are presented in Fig. 03] 
(top panel). The constant slope of Ak b indicates that 
the dust is uniformly distributed between about 2 and 
9 kpc. Because there are no data for D < 2 kpc, these 
interstellar reddening data cannot support or refute the 
claim for the local origin of the dust towards Cen X-3. 
However, by dividing the difference in extinction between 
neighboring bins by the distance between them, one can 
visualize the distribution of extinction elements along the 
sightline (Figure Q31 bottom panel) for D > 2 kpc. The 
results for this region of space are consistent with the 
dust being uniformly distributed, although we point out 
that the error bars allow for a slight increase in dust den- 
sity for D > 3.5 kpc, which is coincident with a slight 
increase in the upper limits to the size of the dust clouds 
shown in Fig. [5] for x > 0.65. 

5. DISCUSSION 

In this paper, we have discovered that most of the dust 
along the line of sight to Cen X-3 is located relatively 
nearby the Sun, and that the distance to the binary sys- 
tem may be smaller than the commonly assumed distance 
of 8 kpc. Although our distance estimate is not particu- 
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larly robust, we note that the best-fit distance and 8 kpc 
are only marginally inconsistent at the 68% confidence 
level, and the best-fit distance is consistent with 8 kpc 
at the 90% level. 

Despite the uncertainty in the distance to Cen X-3, by 
all indications the dust along the line of sight to the sys- 
tem is heavily concentrated within about 300 pc of the 
Sun. Presumably, this dust is from the local Orion spur 
where the Sun resides. Not only do the X-ray halo mea- 
surements, the CO emission measurements, and the den- 
sity of stars along the line of sight support this idea, but 
the halo light curves during egress clearly show a rapid 
response to changes in the point source flux (Fig. [9]). In 
fact, model light curves for a uniform dust distribution 
(not shown), or for a distribution where the primary dust 
cloud is x = [0.20, 0.25], grossly misrepresent the data. 

It is worth examining the difference between the energy 
dependence of the scattering optical depth predicted by 
the WD01 interstellar grain model, and the empirical 
curve describing Cen X-3 as well as data from Cyg X-l 
and GX 13+1 (Fig. HJ). Clearly, the shape of the WD01 
curve and the Cen X-3 empirical curve differ substan- 
tially. One possible explanation for the discrepancy is 
our choice of Ry =3.1. After all, higher values of Ry 
are observed in dense clouds, and we have found that 
most of the dust along the line of sight is located in a 
relatively small region of interstellar space. Yet this is 
not a viable explanation because the difference in the 
scattering optical depth for different values of Ry (over 
2-5 keV) is more in the normalization and not the shape 
of the curve. For a given optical depth at 2 keV, the in- 
crease in optical depth at 5 keV for Ry = 4.0 (Ry = 5.5) 
is 6-7% (14-15%), whereas a -60-80% difference is re- 
quired to explain the energy dependence of r sca in Cen 
X-3. Furthermore, the energy scaling of the scattering 
optical depth in Cyg X-l and GX 13+1 more closely re- 
sembles the empirical curve for Cen X-3 than the WD01 
model curve. In his comparison, Draine (2003) focused 
on the apparent differences between the predicted and 
measured values of T sca / Ay for various X-ray sources, 
and not on the shape of the energy dependence. Our 



result s, and those of ISmith et al.1 (|2002f) and lYao et all 
( 2003), suggest that the WD01 model fails to accurately 
reflect the energy dependence of the scattering optical 
depth. Relative to the WD01 model, there appears to 
be less scattering at low energies, or alternatively, more 
scattering at high energies (due to uncertainties in the 
absolute normalization of the scattering optical depth) 
in Cyg X-l, GX 13+1, and Cen X-3. 

6. CONCLUSIONS 

1. The vast majority of the dust along the line of sight 
to Cen X-3 is located within 300 pc (x < 0.05) of 
the Sun. The most likely location of the dust is the 
local Orion spur. However, we are not sensitive 
to dust located in the vicinity of Cen X-3 because 
data within —3" of the central source are affected 
by pile-up. 

2. The geometric distance to Cen X-3 is 5.7 ± 1.5 kpc 
(68% confidence level), although this result is not 
valid if a large fraction of the dust is concentrated 
within Ax = 0.01. 

3. The energy scaling of the scattering optical depth 
predicted by the WD01 interstellar grain model 
does not accurately represent the results deter- 
mined by X-ray halo studies of Cen X-3. Relative 
to the WD01 model, there appears to be less scat- 
tering at low energies or more scattering at high 
energies in Cyg X-l, GX 13+1, and Cen X-3. 
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APPENDIX 

Chandra/UETG Point Spread Function 

The excellent angular resolution of the Chandra High Resolution Mirror Assembly (HRMA) and the good spectral 
resolution of the Advanced CCD Imaging Spectrometer (ACIS) provide the best opportunity to date for the study of 
X-ray halos. However, careful subtraction of the Point Spread Function (PSF) is still required to produced accurate 
results. The raw images of X-ray sources are the convolution of the X-ray halo surface brightness distribution and the 
PSF of the mirror assembly, modified by the instrumental response of the CCDs, the reduction of effective area at 
larger off-axis angles due to vignetting, and the dithering of the telescope throughout the observation. Exposure maps 
account for all of these effects, with the exception of the PSF. At small halo angles (< 5") MARX 2 simulations can 
be used to model the PSF, but at larger halo angles the results of MARX simulations appear to be unreliable (Gaetz 
2004). Furthermore, with piled- up sources it becomes difficult to obtain the correct PSF normalization. 

In order to address these issues, we chose to develop a representation of the PSF by utilizing all Chandra/KETG 
observations of PKS 2155-304 and Her X-l. These two sources are both out of the plane of the Galaxy and have 
interstellar hydrogen column densities of Ah < 10 20 cm~ 2 along the line of sight. Presumably, this also means that 
the amount of dust along the line of sight is nearly negligible, meaning the radial profiles of the zeroth-order images 
are almost purely due to the PSF, with only a very small contribution due to dust-scattering. We use the HETG 
observations of these sources for two reasons: (1) The grating allows for flux and spectral measurements, making it 
possible to accurately normalize the PSF as a function of energy; and (2) the use of the grating will provide consistent 
analysis when we apply the resulting PSF parameter values to our observations of Cen X-3 (in case the diffraction of 

2 http://space.mit.edu/ASC/MARX/ 
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TABLE 1 

Chandra/HETG PSF Parametrization 



Energy Range 


Normalization a 














(keV) 


(10~ 5 arcscc~ 2 ) 






r 






xl 


1.0—1.2 


9.9 ± 0.3 (3.67o) 


2 


62 


± 





04 


1.40 


1.2—1.4 


9.1 ± 0.4 (4.8%) 


2 


51 


± 





05 


0.99 


1.4—1.6 


7.8 ± 0.5 (6.4%) 


2 


48 


± 





07 


1.24 


1.6-1.8 


7.9 ± 0.6 (8.47o) 


2 


30 


± 





09 


1.92 


1.8-2.0 


8.7 ± 0.5 (5.97o) 


2 


48 


± 





07 


1.21 


2.0—2.2 


10.4 ± U.D (o.07o) 


2 


52 


± 





07 


1.44 


2.2—2.4 


9.U ± (J. 5 (0.870) 


2 


44 


± 





07 


1.20 


2.4—2.6 


7 *7 _i_ n 7 /in /i 0/ \ 
/ . / ± U. f (1U.4/0 j 


2 


36 


± 





12 


2.20 


2.6-2.8 


7.8 ± 0.6 (8.9%) 


2 


12 


± 





09 


1.60 


O Q Q A 

z.o— O.U 


o. / ± U. f (o.o/o ,) 


n 
Z 


on 


1 

± 


U 


1U 




3.0-3.2 


9.6 ± 0.8 (9.7%) 


2 


24 


± 





10 


0.95 


3.2-3.4 


9.6 ± 0.8 (8.9%) 


2 


37 


± 





10 


1.12 


3.4-3.6 


10.0 ± 0.8 (8.7%) 


2 


13 


± 





09 


1.48 


3.6-3.8 


9.6 ± 0.8 (9.3%) 


2 


03 


± 





09 


1.23 


3.8-4.0 


10.5 ± 0.7 (7.4%) 


1 


97 


± 





07 


0.73 


4.0-4.2 


10.8 ± 0.8 (7.8%) 


2 


09 


± 





07 


1.16 


4.2-4.6 


11.8 ± 0.8 (7.2%) 


2 


00 


± 





07 


0.97 


4.4-4.8 


11.7 ± 0.8 (7.7%) 


2 


04 


± 





07 


1.29 


4.6-4.8 


11.6 ± 0.8 (7.3%) 


1 


92 


± 





06 


1.07 


4.8-5.0 


12.1 ± 0.8 (7.4%) 


1 


98 


± 





06 


1.83 



Note. — PSF approximated by co-adding HETG observa- 
tions of Her X-l (ObsIDs 2749, 3821, 3822, 6149, & 6150) and 
PKS 2155-304 (ObsIDs 1014, 1705, 3167, 3706, 3708, 5173, 
6926, 7291, 8380, & 8436). The PSF can be approximated by 
the function S PSF (E,8) = A 5 (E) (0/5")~ r(E \ where A 5 (E) 
is the normalization d" from the point source. Parameters 
apply from 3" to 100" from the central point source; an expo- 
nential cut-off is required for analysis of larger off-axis angles 
(Gaetz 2004). 

a A$(E): Fraction of the point source flux comprising the PSF 
5 /; from the central source in units of 10 — 5 arcsec -2 . The fit 
parameters A$(E) and T are correlated, so the uncertainty in 
A^(E) does not precisely correspond to the uncertainty in the 
integrated PSF flux (/ 2v9S PSF (E , 6)d6) . The percentage in 
parenthesis more accurately represents the fractional uncer- 
tainty in the integrated PSF flux. 

photons from the HETG modifies the PSF). The drawback to this approach, however, is that grating reduces zeroth- 
order effective area and the statistical quality of the images. To try to overcome the problem of limited statistics, we 
co- added the zeroth-order images and first-order MEG/HEG spectra from five observations of Her X-l (ObsIDs 2749, 
3821, 3822, 6149, & 6150) and ten observations of PKS 2155-304 (ObsIDs 1014, 1705, 3167, 3706, 3708, 5173, 6926, 
7291, 8380, & 8436), for a total exposure time of 373 ks. 

For each observation, twenty different radial profiles were extracted in 200 eV energy bands from 1-5 keV using 40 
logarithmically-spaced bins from 0'.'492 (1 pixel) to 110", with spatial filters applied to the readout transfer streak and 
the MEG and HEG diffraction pattern (offset by -4?725 and 5?235 from a line perpendicular to the readout) using 
rectangular regions 10 pixels wide. The radial profiles from each observation were then summed together and scaled 
to the total exposure time, resulting in images with units of [counts/s/pixel]. 

To account for differences in the detection efficiencies between the individual observations, separate instrument maps 
were created for each observation and energy band. Traditional instrument maps account for varying effective area 
and photon detection efficiencies across the CCD chip, but for PSF analysis it is only appropriate to use the effective 
area at the source position (Gaetz 2004). The aspect histograms of each observation were then used to project the 
instrument maps onto the plane of the sky, resulting in effective area maps accounting for bad rows on the CCD and 
the dither of the telescope. The dithered instrument maps, which showed deviations in the effective area at the source 
position of a few percent from observation to observation, were averaged by weighting each individual map by the 
number of counts in each energy and angular bin from the corresponding observation, correcting for the area of the 
spatial filters that were applied to the images. Radial profiles of the exposures maps were then created for the same 
annular regions as the profiles of the zeroth-order images, resulting in maps with units of [cm 2 x counts/photons]. 
The final exposure-corrected radial profiles were then created by dividing the summed counts images by the averaged 
exposure map in each energy band. Using the summed spectrum, fluxes in each energy band were measured in units of 
[photons/cm 2 /s]. By dividing the exposure-corrected images by the fluxes in each energy band, we obtained normalized 
PSF profiles in 200 eV bands. 

To each PSF profile (with background subtracted), we fit the function 

S PSF (E,0)=A 5 (E)^j , (1) 
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where A§ [E) is the normalization 5" from the central point source in units of flux fraction per square arcsecond. When 
modelling the PSF at larger off-axis angles (e.g., 8 > 100") it is necessary to include an exponential cut-off in the 
fitting function to account for the diffuse PSF wings (Gaetz 2004), but at smaller angles (9 < 100") a simple power- law 
representation appears to be adequate. The final product that we obtained through this analysis is the energy- and 
angle-dependent fraction of the source flux comprising the PSF, as shown in Table 1. These results can be safely 
applied to Chandra/HETG observations 3" to 100" from the central source. 
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